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Abstract 

A mesoscopic model of a diblock copolymer is used to study the motion of 
a grain boundary separating two regions of perfectly ordered lamellar struc- 
tures under an oscillatory but uniform shear flow. The case considered is a 
grain boundary separating lamellae along the so called parallel orientation (with 
wavevector parallel to the velocity gradient direction) and along the transverse 
orientation (wavevector parallel to the shear direction). In the model consid- 
ered lamellae in the parallel orientation are marginal with respect to the shear, 
whereas transverse lamellae are uniformly compressed instead. A multiple scale 
expansion valid in the weak segregation regime and for low shear frequencies 
leads to an envelope equation for the grain boundary. This equation shows 
that the grain boundary moves by the action of the shear, with a velocity that 
has two contributions. The first contribution, which arises from the change in 
orientation of the grain boundary as it is being sheared, describes variations 
in monomer density in the boundary region. This contribution is periodic in 
time with twice the frequency of the shear and, by itself, leads to no net motion 
of the boundary. The second contribution arises from a free energy increase 
within the bulk transverse lamellae under shear. It is also periodic with twice 
the frequency of the shear, but it does not average to zero, leading to a net re- 
duction in the extent of the region occupied by the transverse lamellae. We find 
in the limit studied that the velocity of the grain boundary can be expressed 
as the product of a mobility coefficient times a driving force, the latter given 
by the excess energy stored in the transverse phase being sheared. 



1 Introduction 



Self-assembly of block copolymers is one possible route to the development of nanos- 
tructured materials, either directly or as templates. The major challenge that needs 
to be overcome for widespread application of these materials is the development of 
long ranged order in the polymer over scales much larger than the wavelength of the 
mesophase. The purpose of this paper is to investigate the motion of a grain bound- 
ary, an extended topological defect separating two lamellar phases that exhibit long 
ranged order but along different directions. Our study is an extension of ref. [T| in 
that we explicitly consider here that the block copolymer is under an externally im- 
posed, oscillatory and uniform shear flow. The study is motivated by the widespread 
use of both steady and oscillatory shear flows to induce long ranged order in lamellar 
phases. 

The system under consideration is a symmetric diblock copolymer slightly be- 
low its order-disorder transition temperature Tqdt (weak segregation regime). The 
equilibrium phase is a lamellar structure in which nanometer sized layers rich in A 
or B monomers alternate in space. When the copolymer is quenched from a high 
temperature to a temperature T < Tqdt, a transient polycrystalline sample results 
comprised of an ensemble of locally ordered grains but of arbitrary orientations. A 
large number of defects are present in the sample in addition to grain boundaries that 
include dislocations and disclinations. 

Different methods of sample alignment are being investigated experimentally, in- 
cluding substrate induced patterning step induced orientation of thin films 

electric fields that take advantage of a non uniform dielectric constant OIZI, or of 
the existence of ions in the copolymer jB] , and oscillatory shear flows in bulk samples 
Uni HH 1121 Uni mi , We focus here on the latter case as there is no agreement at 
present on the issue of orientation selection as a function of the physical properties 
of the copolymer and the parameters of the flow. For the purposes of the discussion, 
three basic orientations of the lamellae relative to the shear flow are conventionally 
defined: parallel, in which the lamellar planes are parallel to the flow velocity; trans- 
verse in which the lamellar normal is parallel to the flow, and perpendicular in which 
the lamellar normal is parallel to the vorticity of the imposed flow. A review of 
current experimental phenomenology can be found in ^] • • 

In an attempt to clarify the existing phenomenology about orientation selection of 
lamellar phases, we recently undertook a theoretical stability analysis of such phases 
under oscillatory shear flows with and without viscosity contrast between the mi- 
crophases ^3 ^j. At fixed temperature, we found that there is a finite range of 
shear amplitudes within which periodic lamellar structures along a given direction 
exist. For amplitudes larger than a certain critical amplitude, the lamellar phase 
"melts" into a disordered phase, possibly reforming with a different orientation rel- 
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ative to the shear, the new orientation being within the band of allowed solutions. 
Lamellar configurations within the band of allowed solutions can in turn become 
unstable against long wavelength perturbations. The regions of occurrence of this 
secondary instability were given in refs. fTJ QB] as a function of the orientation of 
the lamellae and the shear rate. Generally speaking, it was found that the region of 
stability of perpendicular lamellae is larger than that of parallel lamellae, and both 
considerably larger than the region of stability of the transverse orientation. Our 
results also showed that the critical mode of instability is typically along the perpen- 
dicular orientation, so that the decay of an unstable region of parallel or transverse 
lamellae would lead, at least initially, to lamellae predominantly oriented along the 
perpendicular direction. These results were interpreted through a geometric descrip- 
tion of the lamellar distortion, suggesting that the emerging mode of instability is 
the one that causes the largest decrease in lamellar wavelength under shear. Fi- 
nally, the results were shown to be fairly insensitive to viscosity contrast between the 
microphases. 

While the results just summarized narrow the range over which particular orienta- 
tions can be in principle observed experimentally, they do not provide an orientation 
selection mechanism among competing, linearly stable stationary states. We there- 
fore turn our attention to dynamical aspects of the competition between coexisting 
orientations in a macroscopic sample, and to orientation selection mechanisms of dy- 
namical nature. We focus in this paper on the motion of a front or boundary that 
separates two ordered regions of parallel and transverse orientations. Parallel lamel- 
lae in the mesoscopic description employed are marginal with respect to the shear, 
and therefore unaffected by the flow. Transverse lamellae, on the other hand, are 
compressed by the shear, a fact that is shown to induce boundary motion toward the 
region of transverse orientation. Therefore parallel lamellae are expected to become 
prevalent over transverse lamellae, even in those ranges of parameters of the polymer 
and of the flow in which transverse lamellae are linearly stable. 

The role that topological defect motion plays on structure coarsening under shear 
flows has already been investigated numerically in ^3 1201 • Either a density functional 
description ([H]), or a cell dynamical system model (1201) have been used to study 
domain coarsening of an initially macroscopically disordered conflguration. For steady 
shears, Zvelindovsky et al. ^H] show that the shear is very effective in speeding up 
the formation of lamellar domains, and argue that the perpendicular orientation is 
most stable. Ren et al. [20] more speciflcally focused on topological defect motion 
and annihilation, and the amplitudes of the oscillatory shear necessary to eliminate 
them to form a perfectly ordered lamellar structure. 
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2 Mesoscopic model equation of a lamellar phase 
under shear 



At a mesoscopic level, and for time scales that are long compared with the relaxation 
time of the polymer chain, a block copolymer melt is described by an order param- 
eter ip{r) which represents the local density difference between the two monomers 
constituting the copolymer. The corresponding free energy was derived by Leibler 
in the weak segregation limit (close to Tqdt) EI], and later extended by Ohta and 
Kawasaki to the strong segregation regime If the temporal evolution of ip occurs 
through advection by a flow field as well as through local dissipation driven by free 
energy reduction, ip obeys a time dependent Ginzburg-Landau equation that in the 
symmetric case of equal volume fraction of the two monomers is given by j23j . 

^ + V ■ VV^ = (-^ + - V^^/-) - 5^. (1) 

All quantities have been made dimensionless, including the advection velocity v, 
and the long range polymer interaction coefficient B. The order- disorder transition 
between a disordered phase {ip = 0) and a lamellar phase {ip ^ 0) takes place at 
Bq = 1/4. For B > Bq, -0 is a periodic function of wavenumber go = l/"\/2- 

The physical system under consideration here is a layer of block copolymer, un- 
bounded in the x and y directions, and being uniformly sheared along the z direction 
(Fig. Q). The layer is confined between the stationary z = plane, and the plane z = d 
which is uniformly displaced parallel to itself with a velocity Vpianc = 'jduj cos(co't)x, 
where x is the unit vector in the x direction. 

We first briefly summarize the results of refs. fr7\ ITS] concerning stationary lamel- 
lar solutions in shear flow. In the weak segregation limit e = {B — Bq)/2Bq <^ 1, the 
solution for the monomer composition can be obtained perturbatively in e, 

V^(r) = 2A{t) cos(q ■ r) + Ai{t) cos(3q ■ r) + . . . (2) 

where A{t) ~ (9(e^/^), and Ai{t) and higher order mode amplitudes are of higher order 
in e. We have defined r = xix-|-a;2y + X3(7sin(ti;t)x + z), a vector that is expressed in 
a non orthogonal basis set which follows the imposed shear, and q = (gi, 52,93), the 
wavevector in the corresponding reciprocal space basis set {gi = x — 7 sin(a;t)z, g2 = 
y, g3 = z}. Note that in this new coordinate system, the wavevector of a perfectly 
ordered configuration is stationary. Three orientations relative to the shear can be 
defined as follows: gs 7^ 0, gi = ^2 = is a purely parallel orientation, q2 7^ 0, gi = gs = 
is a perpendicular orientation, and gi 7^ 0, g2 = gs = is a transverse orientation. 

For constant viscosity and if we neglect flow induced by the lamellae themselves, 
the velocity field in the layer is independent of monomer composition, and is given 
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(3) 



dA 
Itt 



a[q\t)]A-M\t)A\ 



(4) 



with q^{t) = q\ + Yl sm{ijjt)qi — ^3]^ + g| and cr{q^) = q^ — q'^ — B. This equation 
can be integrated to give the marginal stabihty boundaries, and the function A{t) 
itself ^7]. From this analysis, a critical strain amplitude 7c was identified, function 
of the orientation q but independent of the frequency cj, such that for 7 < 7c the 
uniform lamellar structure oscillates with the imposed shear, but for 7 > 7c A{t) 
decays to zero; i.e., the lamellar structure melts, according to the terminology used 
by experimentalists. 

The stability of this base lamellar pattern was then addressed by Floquet analysis. 
Regions of stability were obtained for lamellar solutions of arbitrary orientation, that 
were generally largest for orientations near the perpendicular direction, and smallest 
in the vicinity of the transverse direction. As discussed in the introduction, this 
stability analysis provides some guidance on the issue of orientation selection, but we 
wish to extend here the analysis of existence and stability to possible selection by 
dynamical mechanisms. The specific case considered in this paper is the motion of 
a grain boundary separating regions of uniform parallel and transverse orientations 
under oscillatory shear. 

We use in what follows a different form of the equation governing the evolution of 
the monomer composition known as the Brazovski equation (or Swift-Hohenberg 
equation in the fluids literature) 123 121] • Both this equation and eq (H)) lead 
to the same amplitude or envelope equations near onset [27| |2H|, and hence lead to 
identical results in the limit addressed in this paper. The Swift-Hohenberg equation 
for a scalar order parameter is. 



where all quantities are dimensionless. In this units go = I5 although we will retain 
the symbol go iii the equations that follow for clarity of presentation. 

As was the case in the analyses presented in refs. |17| ITH] . we introduce a new 
frame of reference in which the velocity vanishes. In the case of an imposed oscillatory 
shear of amplitude 7 and angular frequency tu, we define a set of non-orthogonal 
coordinates x' = x — a{t)z and z' = z where a{t) = 'jsin{ujt). We assume that 
the system is uniform in the third direction, and therefore simply focus on a two 
dimensional case. Equation (0) transforms to. 
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e^-(V'2 + g2)V-^', 



(6) 
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where 

V" = (l+aW')^-2aW^ + |,. 

A solution of the hnearization of eq ® can be found by assuming 

V'(r') = A cos (q' ■ r') . 
Given that (V + go) = (-g(^)^ + ^o) ^ ^i^h q{tf = ql + (a(t)g^/ - q^,f, we find 

^-^ = eA-{-q{tf + qlfA = a{t)A- (7) 
The disordered solution A = becomes unstable when 



a{t')dt' > 0. 

In analogy with the case analyzed in ref. ^7], we find several instability modes and 
associated thresholds, 

• ga; = 0, ^2 = go e = parallel mode 

• g^ = 0, g^. = ^J ^Jllt^^^Q go e = 3-^4+8^^':^ +8 transverse mode 

• qx = JjE^ ^0, Qz = Ji$^ go e = mixed mode 



3 Envelope equation for a grain boundary under 
weak shear 

We focus on a special configuration comprising two perfectly ordered lamellar do- 
mains, initially oriented perpendicular to each other that meet at a grain boundary. 
We assume that both domains are initially of same wavenumber go at least far from 
the boundary, so that a planar grain boundary would be stationary in the absence of 
shear. We will neglect in this study any back flow induced by the lamellae themselves 
(through osmotic stresses), so that the velocity field v in eq (0) equals the imposed 
shear flow. A schematic representation of the configuration under study is shown in 
Fig. 121 We denote by B the lamellae that lie parallel to the flow field, and note that 
the order parameter ip in this region is unaffected by the flow. Transverse lamellae 
are denoted by A. If the A lamellae were to adiabatically follow the imposed flow, 
both orientation and wavelength would be a periodic function of time as illustrated 
schematically on Fig. |21 Because a local change in wavenumber away from go always 
leads to a free energy increase in region A, while the free energy in region B remains 
unchanged, we anticipate grain boundary motion from region B to region A, thereby 
increasing the area occupied by parallel lamellae. 
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We derive a set of amplitude equations from eq (jH)) by using a multiple scale 
approach. For e <^ 1, it is possible to extract the slow evolution of both the lamellae 
and the grain boundary (on a slow time scale et) by expanding the order parameter ip 
in both regions around a periodic function, with amplitudes that are slowly varying in 
the grain boundary region (of very large extent in this limit). Our derivation follows 
closely that of Tesauro and Cross for the case of no flow 29]. 

The analysis is restricted to shears of small amplitude and low frequency. Specif- 
ically, since —q{tY + = —ci^ql' + '^^aqx'qz', the shear contributes to eq ((Tj) with two 
terms, one of order 0{a^) and the other of order 0{a^). Consistency with the expan- 
sion in the weak segregation limit requires that {—q(tY + q^) ~ C(e), a requirement 
that dictates the magnitude of the shear amplitude. For the case considered below 
involving a grain boundary between parallel and transverse orientations, the cross 
derivative term —2aqx'qz' acting on the reference state vanishes, and therefore we will 
have (-g(t)2 + g2)2 _ ^(^4) ^ _ 0{e^/''). 

We start by assuming that the slowly varying amplitude of mode e*''"^' in region 
A has as characteristic length scales X = e^^^x' and Z = e^^^z', while the mode 
g«(?o2 domain B has characteristic scales X = e^^^x' and Z = e^l'^z' . We further 
assume a weak shear so that a ~ e^l^. This scaling is appropriate for an initial 
configuration in which the transverse lamellae are parallel to the grain boundary (and 
therefore in the limit of small shear amplitude, the transverse lamellae will remain 
almost perpendicular to the parallel lamellae for all times). If the angle between the 
transverse lamellae and the grain boundary is finite, then the scaling a ~ e^/^ needs 
to be introduced instead. The operator (V'^ + gg)^ can now be expanded in powers 
of e as 

+e=^/^(2LoL3 + 2L1L2) + e(2LoL4 + 2L1L3 + L^), (8) 

where we have defined, 

Lo = dl + dl + ql 

Li = 2{dx'dx + d^idz - ad^'d^') 

L2 = d\ + dz + '^{dx'dx + d^idz - ad^/dx - ad^>dz) + a^d^, 

L3 = 2[dxdx + dzdz - a{d^>dx + d^>dz + dxdz) + a^d^tdx] and 

U = dl + dl- 2a{dxdz + dxdz) + a\2d,,dx + 4). 

We also expand as 

^ = ei/Vo + e^i + e'/'^2 + - (9) 
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and assume that both the frequency of the imposed shear, and the associated variation 
of ip is over a slow time scale T = et. From eqs (jHI), (0) and Q we obtain, at 0(e^/^) 
the equation 

- Ll^P^ = 0, (10) 

which admits the solution 

^0 = ^[Aoe*'?""' + Boe'««"' + c.c], (11) 
v3 

with Aq and Bq functions of X, X, Z and T. At 0(e), eq (jHl) reduces to 

- Llij^ = Ll^o, (12) 

where we have used the fact that Lq^/'o = 0, and taken advantage that the cross 
derivative term ad^'dz' vanishes when acting on ipo, so that the solution at this order 
is also time independent (in the sheared frame of reference). Since ipo is an eigenmode 
of Lq with zero eigenvalue, the right-hand side of eq (fT^ must vanish in order for 
it to admit a solution. Solvability requires that the scalar product (■?/'(]' iL^-^o) = 0, 
that is, the right hand side of eq (fT^ is orthogonal to the the zero eigenf unctions 
of the adjoint of Lq. But (ipolLlipo) = (LfipQlLixpo) = ULi-^oll^ = 0, from which 
it follows Liipo = 0. As a result of this solvability condition, Aq and Bq must be 
independent of X and Z, respectively. Equation (fT2jl then reduces to eq (finjl and 
V'l = [Aie^-^o^' + 5ie^««"' + c.c,]/v^. 

Finally, at 0(e'^/^), the multiple scale analysis yields the equation 
Llip2 = -dri^o + i^o- i^o - (2^o^4 + 2L1L3 + LD^Jq - (2L0L2 + L?)^i. (13) 

Again, the functions ipo and ipi are zero eigenmodes of the operator Lq, so that the 
projections of the terms in the right-hand side of eq (fT^ on these eigenfunctions must 
vanish. From this condition, we obtain the following amplitude equations 

drAo = {1 - [2tqo{dx - adz) + 9| - qyf}Ao - \AofAo - 2\BofAo, (14) 

and, 

OtBo = {l- [2tqoidz - adx) + 4]'}^o - l^o^fio - 2\Ao\^Bo, (15) 

where we have used the fact that Lqi/jq = Liipo = and have set Llipi = 0. This 
set of equations (fT^ and (fT3j) governs the evolution of the slowly varying envelopes 
of the base lamellar pattern, including variations both in the direction parallel and 
perpendicular to the grain boundary. We now restrict our attention to the case of a 
planar grain boundary; hence we do not consider any dependence of the amphtudes 
on the coordinate parallel to the grain boundary. Transverse perturbations of the 
grain boundary are expected to decay back to planarity, and such decay will not 
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be considered here. Therefore Aq and Bq depend only on X and X respectively. 
Equations (fT^ and (fT^ simplify to 

OtAo = [1 - {2tqodx - qyf]Ao - |Ao|Mo - 2|5o|Mo, (16) 

and, 

OtBo = [1 - (4 - 2iaqodxY]Bo - |5op5o - 2|Aop5o. (17) 

In the absence of shear {a = 0,x' = x), these two equations reduce to the case studied 
by Manneville and Pomeau jSHI, and by Tesauro and Cross 

Finally, we let c^^'^Aq = r^e^'^^, g^^^Bq = rse**^-^ where to lowest order 0^ and 0^ 
are independent of X and X j^, and re- write the set of amplitude equations in the 
original (unsealed) set of variables. The resulting equations for and read 

dtTA = ^qldlvA + (e - qtayA -r\- 2r|rA, (18) 

and, 

dtTB = -dliTB + ^a^qldliTB + trs - r% - 2r\rB- (19) 

This is one of the main results of our calculation. At this order, there are two 
contributions to the amplitude equation arising from the shear. One in eq ()19|) which 
multiplies the second normal derivative of the amplitude r^. This contribution is 
nonzero only in the grain boundary region, of twice the frequency of the imposed 
shear, and describes variations of the amplitude of the parallel lamellae due to the 
changing orientation of the grain boundary with respect to the lamellar planes. As we 
will show below, this term alone originates an oscillatory contribution to the velocity 
of the grain boundary of zero average. The second contribution is the term qod^TA in 
eq (jl8j) . This term does not contain a spatial derivative, and therefore is important 
in the entire bulk region A where the amplitude does not vanish. It leads to a 
change in the amplitude of the uniform transverse lamellae as they are advected by 
the flow. The corresponding change in the free energy of region A does not average 
to zero over a period and is responsible for the net motion of the grain boundary, as 
shown below. 

The amplitude equations (fTS|) and (fT^ need to be supplemented with appropriate 
boundary conditions. First we have that r^(— oo,t) = and rB{+oo,t) = 0. Fur- 
thermore, at large distances from the grain boundary inside domain B, reduces to 
the constant ^/e, independent of the flow parameters. By contrast, the amplitude 
inside domain A satisfies the equation dtr^ = {e — qQa'^)rA — r\ in the limit x' — > -|-oo. 
That equation admits the solution 

r^(+oo,t) = {e-^W[2 f e^^^'Ut' + rf{+oo,0)]y'/^ (20) 

Jo 

where f{t) = (2e - |g^7^) t + g^7^ {^-^^ - ^ij^). The asymptotic behavior of eq 
at large times changes qualitatively with the sign of the constant 2e — |(?o7^- 
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When this constant is negative, the prefactor e"-'^^*-' diverges exponentially with time 
and rA(+oo, t) decays to zero. If, on the other hand, 2e — |gQ7^ > 0, Laplace's method 
can be used to approximate the integral in the expression for r^, which reduces to a 
periodic function 



rA{+oo,t) 



9{t) 



h{t) 



-1/2 



(21) 



where g{t) = a;gQ7^(sin 2ujt — \ sin Aujt) and h{t) = 2e + gQ7'^(cos 2u}t — \ cos Aut — |). 
The condition 2e — |go7^ = which separates these two cases can be understood 
in terms of a maximum strain amplitude 7* = (Se/Sgg)^''^ above which the lamellar 
phase of domain A will melt. Note that 7* is independent of the shear frequency uj. 



We have numerically solved the coupled, one dimensional equations (jl8|) and (jl9|) . 
The results that will be shown correspond to e = 0.04, = 1 and 7 = 0.3, and a 
variety of shear frequencies. In the calculations, region A was surrounded by two 
identical domains of parallel (B) lamellae so that periodic boundary conditions in 
both directions could be used. The equations were integrated with a pseudo-spectral 
algorithm described in detail in ref. ^Tj. We used a computational domain of size 
L = 4096 and a grid spacing As = 0.5. The time interval At was chosen a function 
of the period of the shear as At = 277/(5000007). Stationary solutions obtained in the 
absence of shear provided the initial conditions for and tb- 

The instantaneous location of the grain boundary x* was defined as rsix*) = 
\/e/2, and its velocity Vgb as the rate of change of x*. FigureElshows 

Vgb cL function of 

time for several values of the angular frequency u. Time has been scaled by the period 
r = 271 /u of the applied shear. Positive (resp. negative) values of Vgb indicate motion 
toward lamellae in domain A (resp. B). Following an initial transient, the velocity 
oscillates in time at half the period of the shear, with an amplitude that decreases with 
the frequency. Note also that the average velocity is positive; i.e., motion is directed 
toward domain A. It is possible to further interpret these results by obtaining an 
analytic approximation for the boundary velocity in the limit of very low frequencies. 
If the frequency is sufficiently low, the order parameter ip (or the amplitudes and 
Tb) will adiabatically follow the motion of the boundary. In this limit it is possible to 
invoke a quasi-static approximation according to which rA{x',t) ~ r\{x' — x'g^{t);a) 
and tb^x', t) ~ r^(x' — x'^f^it); a), where r\ and are stationary solutions of eqs (fT5|) 
and (fT^ (with the boundary conditions given, including eq (EH) that still formally 
depend on the parameter a. This latter dependence results from the dependence of 
the stationary amplitude profiles r\ and on the state of shear of the system given 
by the parameter a. Then {dtrA^dtrs) — —Vgb{dxrA,dxr%) where we have introduced 
the notation x = x' — x^^. Following Manneville and Pomeau jHD], we multiply eq ()18|1 
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by dxr\ and eq (fT^ by ^^r^, add the two equations and integrate over x to obtain 

/ + 00 
dx [{d^r\f + [d^r^B?] = K{+oo) - K{-oo), (22) 
-oo 



where 



+4a'ql{d,r%f - 2 



The integral in the left hand side of eq ()22|1 is an inverse mobility or friction coefficient, 
whereas the right hand side is the effective driving force for grain boundary motion. It 
is equal to the static free energy increase of the configuration upon shearing relative to 
the planar, unshared boundary and can be evaluated by using the values r\{—oo) = 
0, r^(— oo) = ^/e and r^^l+oo) = 0. Furthermore in the quasistatic limit the function 
g{t) appearing in eq (PT|) is small, and r% ^ ^/7^^q^. As a result, K{—oo) = 
6^4:, K{+oo) = (e - qyY/4: and 

^ ^ go7^ sin^(a;t) [2e - q^-f^ sin^(^t)] 
4 J^^ dx mr\y + {d,r%y] ■ 

We now compare this result with those obtained by numerical integration of the 
governing equations and shown in Fig. El for a range of angular frequencies. In order 
to compute the inverse mobility coefficient that appears in the denominator of eq (j23|) 
we have obtained rA{x',t) and rB{x',t) directly by integration of eqs (fTH|l and (fTUj) . 
Equation (PH)) agrees very well with the numerical value of Vgb at very small shear 
frequencies, as illustrated in Figure El for u = 0.001. The agreement progressively 
deteriorates as the angular frequency increases. 

We finally note that although the time dependence of Vgb changes significantly with 
uj, its average over a period (vgb) = (l/r) v{t) dt only shows a weak dependence on 
shear frequency, as shown in figured We find, for example, that increasing the shear 
frequency by two orders of magnitude causes a decrease of only 5% in the average 
speed of the grain boundary. This result follows from the fact that the effective 
driving force K{—oo) — K{+oo) responsible for the motion of the grain boundary 
is independent of the shear frequency in the quasistatic limit and increases only 
marginally at larger frequencies. As a result, variations in velocity with uj arise solely 
from the small changes of the inverse mobility coefficient on frequency. Somewhat 
unexpectedly, the quasistatic approximation does quite well at predicting the average 
velocity of the grain boundary for a wide range of frequencies. We show in Fig. HI 
(vgb) obtained by averaging both sides of eq over a period, and the numerical 
results of Fig. |21 also averaged over a period of the shear. 
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4 Discussion 



Following a quench of the diblock copolymer from an initially disordered configura- 
tion at T > ToDT to a final temperature below Tqdt, the following qualitative picture 
emerges concerning the asymptotic, long time selection of a lamellar orientation rela- 
tive to the shear. In the absence of shear (7 = 0), initial composition fluctuations are 
amplified exponentially, with a growth rate that is isotropic. Lamellar regions emerge 
and coarsen as a function of time. Coarsening rates and the role of topological de- 
fects in a two dimensional system have been discussed in ref. [32] • To our knowledge, 
a similar investigation in three dimensions has not been carried out. If the quench 
takes place under shear, the mean field instability threshold depends on orientation, 
as shown in ref. ^7] . The first threshold is to a mixed parallel-perpendicular mode at 
Be = 1/4, followed by a bifurcation to a transverse mode at Be = 1/4 — 7^/32 + 0(7^), 
and to a parallel-transverse mode at Be = 1/4: — 7^/8 + 0{'y'^). Therefore for shallow 
quenches fluctuations along different orientations would be amplified at different rates 
leading to predominantly parallel and perpendicular oriented domains even from an 
isotropic initial condition. Thermal fluctuations, on the other hand, are known to 
significantly modify these conclusions ^31 1331 1^ l^4j . In particular, thermal fluctua- 
tions render the mean field supercritical bifurcation a weakly subcritical bifurcation, 
with a transition temperature that increases with 7. In any event, the distribution of 
orientations following a quench in shear flow is expected not to be isotropic. 

Regardless of whether the copolymer is quenched in shear flow or not, a macroscop- 
ically disordered configuration will result at intermediate times comprising regions of 
well saturated monomer composition, but with a wide distribution of lamellar orien- 
tations. The distribution of observable orientations is reduced by the shear, as those 
orientations that are unstable against long wavelength perturbations will quickly de- 
cay when the monomer composition locally reaches nonlinear saturation. Insofar the 
melt is Newtonian at the low shear frequencies investigated in refs. [T71ITH] . one would 
expect lamellae that are predominantly perpendicular, and to a lesser degree parallel, 
with small projections of the lamellar wavevector on the transverse direction. There 
also exists a small region of stable transverse lamellae. Further structure coarsening 
under shear will involve an initially anisotropic distribution of orientations, and hence 
is expected to be qualitatively different from the isotropic case. In metals, for exam- 
ple, coarsening of an initially anisotropic distribution of orientations leads to texture 

m 

The results of this paper further indicate that regions where parallel and trans- 
verse lamellae meet will move, even when both parallel and transverse lamellae are 
linearly stable. The net motion of the grain boundary is driven by free energy re- 
duction because parallel lamellae are unaffected by the shear, whereas transverse 
lamellae are elastically compressed, a compression that leads to an increase in en- 
ergy that is relieved through boundary motion. Since a similar argument can be 
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made for a boundary separating perpendicular and transverse lamellae, we would ex- 
pect that regions of a macroscopic sample oriented along any combination of parallel 
and perpendicular orientations will grow at the expense of any remaining transverse 
lamellae. 

We are then confronted with an interesting question concerning the behavior of 
boundaries separating parallel and perpendicular lamellae. Both orientations show 
fluid like response to the shear, in contrast with the (one dimensional) elastic response 
of the transverse orientation. In fact, the flow does not couple to the monomer 
composition for the simple model of Newtonian melt with constant viscosity adopted 
here, and hence there are no shear flow effects on this type of boundary. 

If viscous or viscoelastic contrast between the microphases is allowed, a secondary 
flow appears which is orientation dependent |2ni ^] • The velocity field of this sec- 
ondary flow is parallel to the lamellar planes (assuming incompressibility) , and largest 
for a uniform parallel configuration, while it vanishes for a uniform perpendicular 
configuration. This flow is weak in the weak segregation limit, and has negligible 
consequences on the stability of a lamellar configuration against long wavelength 
perturbations ^H|- However, its possible effect on boundary motion has not been 
investigated yet. For long wavelength modulations of the type described by ampli- 
tude equations, it is possible that the effect of these secondary flows can be subsumed 
into an effective constitutive relation for the dissipative part of the stress tensor as a 
function of the velocity gradient tensor. This constitutive relation has to be compat- 
ible with the uniaxial symmetry of a lamellar phase, in analogy with other uniaxial 
systems such as nematic or smectic liquid crystals. Additional viscosity coefficients 
would enter the constitutive law, as well as an explicit dependence on q, the slowly 
varying normal to the lamellar planes. In this case, the effective viscosity of a region 
of parallel lamellae is different from that of perpendicular lamellae thus leading to 
secondary flows that are asymmetric with respect to the boundary, and possibly to 
boundary motion. This possibility is currently under investigation. 
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FIGURE CAPTIONS 



Figure 1. Schematic representation of the geometry considered including the shear direc- 
tion, and the three different lamellar orientations discussed in the text. 

Figure 2. Schematic representation of a planar grain boundary that separates regions of 
parallel and transverse lamellae being uniformly sheared. 

Figure 3. Grain boundary velocity as a function of time obtained by numerical solution 
of eqs (fT^ and ()15|1 . Four different angular frequencies are shown: (in order 
of decreasing amplitude) oj = 0.001,0.01,0.05, and 0.1. Also shown is the 
quasistatic approximation of eq calculated at the lowest angular frequency 
u = 0.001. The curve is indistinguishable in the graph from the corresponding 
numerical solution. 

Figure 4. Temporal average of the grain boundary velocity as a function of the angular 
frequency of the shear. The symbols correspond to the time average of the 
numerically obtained velocities shown in Fig. |21 and the solid line is the time 
average of the quasistatic velocity given in eq (j2Sl)- 
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Figure 1: Schematic representation of the geometry considered including the shear 
direction, and the three different lamellar orientations discussed in the text. 
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Figure 2: Schematic representation of a planar grain boundary that separates regions 
of paraUel and transverse lameUae being uniformly sheared. 
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Figure 3: Grain boundary velocity as a function of time obtained by numerical so- 
lution of eqs ()14|) and ()15|) . Four different angular frequencies are shown: (in order 
of decreasing amplitude) u = 0.001, 0.01, 0.05, and 0.1. Also shown is the quasistatic 
approximation of eq calculated at the lowest angular frequency u = 0.001. The 
curve is indistinguishable in the graph from the corresponding numerical solution. 
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Figure 4: Temporal average of the grain boundary velocity as a function of the 
angular frequency of the shear. The symbols correspond to the time average of the 
numerically obtained velocities shown in Fig. El and the solid line is the time average 
of the quasistatic velocity given in eq (j2Sl). 
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